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ABSTRACT 

Context. Galactic and extragalactic relativistic jets are surrounded by rich environments that are full of moving objects, 
such as stars and dense medium inhomogeneities. These objects can enter into the jets and generate shocks and non- 
thermal emission. 

Aims. We characterize the emitting properties of the downstream region of a standing shock formed due to the interaction 
of a relativistic jet with an obstacle. We focus on the case of red giants interacting with an extragalactic jet. 

Methods. We perform relativistic axisymmetric hydrodynamical simulations of a relativistic jet meeting an obstacle of 
very large inertia. The results are interpreted in the framework of a red giant whose dense and slow wind interacts with 
the jet of an active galactic nucleus. Assuming that particles are accelerated in the standing shock generated in the jet 
as it impacts the red giant wind, we compute the non-thermal particle distribution, the Doppler boosting enhancement, 
and the non-thermal luminosity in gamma rays. 

Results. The available non-thermal energy from jet-obstacle interactions is potentially enhanced by a factor of ~ 100 
when accounting for the whole surface of the shock induced by the obstacle, instead of just the obstacle section. The 
observer gamma-ray luminosity, including the effective obstacle size, the flow velocity and Doppler boosting effects, 
can be ~ 300(7j/10) 2 times higher than when the emitting flow is assumed at rest and only the obstacle section is 
considered, where 7 j is the jet Lorentz factor. For a whole population of red giants inside the jet of an active galactic 
nucleus, the predicted persistent gamma-ray luminosities may be potentially detectable for a jet pointing approximately 
to the observer. 

Conclusions. Obstacles interacting with relativistic outflows, for instance clouds and populations of stars for extragalactic 
jets, or stellar wind inhomogeneities in microquasar jets and in winds of pulsars in binaries, should be taken into account 
when investigating the origin of the non-thermal emission from these sources. 
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1. Introduction 


Relativistic jets are surrounded by rich environments, and 
stars and dense medium inhomogeneities can enter into 
these jets leading to th e formation of shocks, a s cons i dered 
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are, 


some cases, almost standing shocks, at rest in the labora¬ 
tory frame (LF). A similar situation can occur in the termi¬ 
nation region of other re lativistic outfl ows, s uch as p ulsar 
winds in binary systems (Paredes-Fortuny et al. 2015). 

In several previous papers that modelled the radiation 
from relativistic shocks at rest in the LF, the emitting re¬ 
gion was circumscribed to the closest obstacle vicinity and, 
as the flow there is only mildly relativistic, the radiation 
Doppler boosting was not considered. However, there could 
be more jet energy available for non-thermal particles than 
assumed in these works, as the effective interaction area 
may be much larger than just the obstacle section. In ad- 
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dition to that, the flow initially moves at ss c/3 right after 
the shock, and then is quickly accelerated further down¬ 
stream because of strong pres sure gradients. This eff ect i s 


seen for instance in fig. 2 in Bosch-Ramon et al. (2012). 


Therefore, even relatively close to the obstacle at rest, the 
shocked flow velocity is high enough for Doppler boosting 
to significantly enhance radiation in the emitting fluid di¬ 
rection of motion. The motion of the shocked flow could be 
also considered for an ob stacle that is itself accelerated by 
the jet ram pressure (e.g. Barkov et al.|2012a Khangulyan 


et ak|2013 1, but this would require accounti ng for the com 

Td 


Bosch-Ramon 


plex dynamical evolution of the obstacle (see 
|et al.|[20l2l for hydrodynamical simulations predicting the 
obstacle disruption). 


In this work, we study the kinetic-to-internal (possibly 
non-thermal) energy conversion, gamma-ray radiation effi¬ 
ciency, and Doppler boosting effects, in the shock and the 
shocked flow of a relativistic hydrodynamic jet interacting 
with an obstacle at rest. The obstacle could be a star or a 
cloud in the case of an extragalactic jet, or a stellar wind in¬ 
homogeneity in the case of a galactic microquasar jet. Since 
this is a relevant and illustrative case, we focus here on the 
case of a population of red giant stars (RG) interacting 
with the jet of an active galactic nucleus (AGN). We first 
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characterize the non-thermal luminosity from the shocked 
jet region accounting for hydrodynamical information, and 
second we adopt a phenomenological description of the spa¬ 
tial distribution of the RG in the regions crossed by the jet 
as it goes through the galaxy. To compute the emission, we 
consider mainly inverse Compton (IC) scattering, assume 
that synchrotron radiation is a minor emission channel, and 
put special emphasis on the GeV-TeV range. Studies such 
as this are needed to characterize the non-thermal emission 
from interactions between stars and jets in AGN, which is 
important if the AGN gamma-ray radiation is to be und er- 
stood (e.g. Barkov et al.||2012a Khangulyan et al.|2013 ). 

Concerning other types of sources, the characterization 
of the radiation properties from jet-obstacle interactions 
could also shed light on the main factors shaping jet vari- 


ability at high energies in microquasars (e.g. Owocki et al. 


2009), or the post-shock emission in pulsar-wind termi¬ 


nation shocks in some binary systems (e.g. Bosch-Ramon 


2013). Regarding jets with strong magnetic fields with po¬ 


larity changes, it has been proposed that the interaction 
with obstacles may lead to mag netic reconnection an d effi¬ 
cient non-thermal emission (see Bosch-Ramon||2012 ). This 
kind of interaction is very complex given the strong dynam- 
ical role of the magnetic field and its non-trivial topological 
properties and deserves a specific treatment, which is be¬ 
yond the scope of this paper. 


2. The jet-RG wind interaction 

2.1. Physical scenario 

The study presented here is based on the interaction of a 
jet with an RG surrounded by its dense and slow wind. The 
basic hydrodynamical ingredients of this study are: a highly 
relativistic, supersonic flow with parallel fluid lines on the 
scales of interest; a cold obstacle at rest; and a fraction of 
the internal energy downstream of the standing shock being 
in the form of non-thermal particles. The initial speed of 
the RG is taken to be <C c if it is to be approximated as at 
rest, so the problem can be considered axisymmetric around 
an axis crossing the star and with the jet flow direction. In 
addition, the RG wind must have enough inertia to stand 
the jet impact while it penetrates into the jet, which can 
be expressed as p rg v^ g > Lj/nR^c, where p Ig and v rg are 
the RG wind density and velocity perpendicular to the jet, 
and Lj and Rj are the jet power and radius, respectively. 
In addition, in the standing shock scenario the velocity ac¬ 
quired by the RG wind due to jet impact must be -C c. If 
these conditions are fulfilled, the problem becomes two di¬ 
mensional, and the steady state of the physical system can 
be taken as the steady state of the shocked jet flow, which 
significantly lightens the hydrodynamical calculations pre¬ 
sented in what follows. 


2.2. Simulations 


The jet-RG wind interaction was simulated in 2D assum¬ 
ing axisymmetry. We solved the equations of relativistic 


hydrodynami c s usi ng a Marquina Riemann solver (Donat 
fe Marquina| 19961, second-order spatial reconstruction 
scheme (minbee slope limiter, plus an additional dissipa¬ 


tion term from Colella 1985), and an adiabatic gas with 
index 7 = 4/3. The resolution is 600 cells in the vertical 
direction, the 2 -axis, and 300 in the radial direction, the 
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Fig. 1. Distribution of density in the flow frame in cgs units. 


r-axis, and the physical size is z max = 2 x 10 15 cm in the 
2 -direction, and r max = 10 15 cm in the r-direction. The jet 
power injected in the grid is ss 10 37 erg s -1 , with a Lorentz 
factor 7 j = 10. We fix the RG wind radius to R Ig = 10 14 cm. 
A snapshot of the density distribution in the steady state is 
presented in Fig. |T] The shocked jet material is clearly seen 
in light blue, forming a cometary tail pointing upwards and 
surrounding the RG wind, red with green borders. 

For simplicity, we modelled the RG wind region as a 
sphere with a density 10 8 times that of the jet in the FF. 
As the shocked RG wind is much slower than the shocked 
jet flow, we considered it static, and thus the steady state of 
the simulation is quickly reached. In a more realistic setup, 
there would be more wind momentum flux perpendicular to 
the jet direction, slightly widening the shocked flow struc¬ 
ture, but this setup would have required a simulation ~ 10 3 
times longer. 

We chose the grid resolution such that numerical vis¬ 
cosity played a negligible role in the generation of internal 
energy in the shocked flow. There is some mixing in the RG 
wind-shocked jet flow boundary, but it has a minor impact 
on the radiation estimates presented below. This mixing 
slows down the shocked jet flow close to the axis, reduc¬ 
ing the effect of Doppler boosting and thus weakening the 
emission from those (small) regions. 

For illustrative purposes, the simulation result is inter¬ 
preted as an AGN jet of Lj « 10 44 erg s _1 and Rj = 1 pc, 
much larger than the grid, interacting with an RG wind re- 
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gion whose size is derived assuming jet-stellar wind pressure 
balance 


Rr 


M w v w c/ALj x i?j ~ 10 i 4 cm, 


( 1 ) 


for a stellar mass loss M w = 10 ~ 8 Mg/yr and a wind ve¬ 
locity v w ~ 2 x 10 7 cm s _1 , typical of an RG. Note that 
for a real RG wind the actual size of the obstacle would 
be slightly larger than for the case of an obstacle of very 
large inertia, as the shocked wind region should be added 
in Eq. ([l]). 

Relativistic jets usually present half semi-opening angles 
of several degrees, cq ~ 0.1 rad, which for a supersonic 
jet without lateral causal co ntact yields 7 j = l/ay ~ 10 
(see Khangulyan et al. 2013 and references therein). This 


condition implies an interaction height of Zj = 7 ji?j ~ 10 pc. 


3. Non-thermal emission 

We use the hydrodynamical information obtained in Sect. [2] 
to derive some important physical properties that allow 
for the characterization of the non-thermal emission. These 
properties are the radiative efficiency, the distribution of 
the non-thermal energy, the Doppler boosting, and the ef¬ 
fective section of the RG wind as an obstacle for the jet. 


3.1. Radiative efficiency and non-thermal energy distribution 

To estimate the non-thermal emission, the radiative effi¬ 
ciency of the shocked flow is needed. This efficiency can be 
quantified as the radiative (l/t ra d) to total (l/i c ) cooling 
rate ratio (/ rad ): 


/rad — AadA. 


^radA^rad+Arad) (1+tradAnrad) > (2) 


where l/t c also takes non-radiative cooling, l/t mad , into 
account. The timescale for non-radiative cooling, i.e. par¬ 
ticle escape and adiabatic losses, is estimated in the flow 
frame (FF) as the distance between each point to the RG 
location (d rg ) divided by the local velocity (v s h) times the 
shocked flow Lorentz factor ( 7 s h): tnrad = d rg /t? s h 7 sh- For 
the radiative timescale we adopt: t rad = Cic u~ 1 E~ 1 , where 
Cic ~ 26 in cgs units. This timescale is derived assuming 
IC in Thomson, and taking a reference electron energy of 
Eic = m^c^/eo, where Pic is roughly the energy of the 
IC cross-section maximum, and eo is the target photon en¬ 
ergy and u the energy density of a photon field coming 
from the RG, all in the FF. An effective stellar temper¬ 
ature of 3000 K has been taken, typical for an RG. The 
star luminosity, needed to obtain it, has been determined 
through L* = M w v w c (LF). This is a rather conservative 
assumption, as in fact only a fraction of all the star radia¬ 
tion momentum goes to the stellar wind. Note that protons 
are not taken into account as emitting particles. 

Figure [2] shows the distribution of / rad in the shocked 
flow. As seen in the figure, for the adopted jet and RG 
wind properties, / rad is everywhere below one for the case 
simulated, i.e. the flow is adiabatic. In this case, the thermal 
pressure (P) of the shocked flow is a good proxy for the 
non-thermal energy distribution, which can be taken to be 
proportional to P all through the shocked flow. In Fig. [3j 
the distribution of P in the shocked flow is shown in a colour 
scale, from red to light green. 
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Fig. 2. Distribution of radiation efficiency for the shocked 
jet region. 


The characterization of the radiative and the non- 
radiative timescales, and of the non-thermal energy distri¬ 
bution, allows for the computation of an approximate total 
luminosity for IC in the FF, Lie, simplifying the IC emis- 
sivity as that at Eic- For moderately magnetized shocked 
flows, the synchrotron luminosity L sync may be higher than 
Lie, as Async oc P^,where a is the power-law index of 
the particle energy distribution, and P max is the maximum 
electron energy, which could be Pic ■ Under efficient par¬ 
ticle acceleration and 7 j > 10 , the synchrotron spectrum 
may contribute to or dominate in the gamma-ray band. 
However, as the magnetic field is not well constrained, we 
will assume here that the synchrotron component can be 
neglected. 

The adopted radiation cooling model is simple, although 
it accounts for the basic IC cross-section properties, and 
the target density behaviour under relativistic motion of 
the emitter. Therefore, the model should yield reasonable 
luminosity predictions for IC gamma rays, if gamma-ray 
pair-creation absorption is negligible. 

3.2. Doppler boosting and obstacle effective radius 

The radiation enhancement due to Doppler boosting is 

= l/7sh(l - cos (9 sh /3 sh ) 4 , (3) 

which can be approximated as « 167 4 h when the angle be¬ 
tween the jet and the observer directions is 0 o bs l/7sh- 

Figure [4] shows the distribution of the factor of enhance¬ 
ment of the emission from the shocked flow because of 
Doppler boosting. As seen in Fig. |ij values of of a few 
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Fig. 3. Distribution of pressure in cgs units for the shocked 
jet region. 
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Fig. 4. Distribution of the Doppler boosting enhancement 
(<5g h ) for the shocked jet region. 


are found right after the shock, but further downstream (5^ h 
achieves values as high as several thousand. 

Another important result from our simulation is derived 
when computing U s h, the rate of the total thermal energy 
generated in the interaction region in the LF, where the 
total thermal energy is the sum of internal energy plus the 
pressure component. The value of U s h has been computed 
as the fraction of Lj that escapes the grid in the form of 
thermal energy. From this, we can derive the effective radius 
of the RG wind when acting as an obstacle for the 
jet 

Csh = ( Res/Ri ) 2 Tj « 4 x 10 36 erg s" 1 , (4) 

which turns out to be R e g ~ 6 R rg because of effective ki¬ 
netic to internal energy conversion even in regions far from 
the RG. This effective conversion manifests itself locally as 
well, as shown by Fig. [5} 

The value of R e g is actually dependent on R rg and 2 max . 
To check this dependence, two additional simulations were 
carried out, both with twice the cell size of the presented 
simulation, one with the same grid size, and the other with a 
grid that is twice as large, and both with the same remain¬ 
ing parameters. These simulations show that i? e ff grows 
slowly with r max - This is expected to happen as the oblique 
shock becomes weaker farther away from the RG, but the 
actual weakening of the shock is sensitive to the specific 
value of the unshocked flow internal energy, which is un¬ 
known. Note also that to derive R c g all the grid is taken 
into account; in reality, the efficiency of energy conversion is 
larger if the grid outermost region without shock presence 


is neglected. Therefore, R e g = 6 R Ig may be considered a 
lower limit. 

4. The 1C luminosity from a jet-obstacle 
interaction 

4.1. The 1C luminosity to the observer 

Using the hydrodynamics results and under the simplified 
emission model described, we can derive an approximate 
observer IC luminosity of the shocked flow within the grid 
from 

r 

Aic.obs ~ / £* h (l - cos6bobs)r?NTTP(7 - l) _1 t“^dn sh , (5) 

J sh 

where dv s h is the shocked flow differential volume, and 0* o bs 
the angle between the stellar photons and the observer di¬ 
rection, both in the FF. The term (1 — cos0* o bs) is a term 
related to the IC scattering probability-dependence on the 
stellar photon and observer directions. The quantity ?7 nt 
is the ratio of non-thermal to total thermal energy den¬ 
sity; it is a free parameter of value between 0 and 1. For a 
typical accelerated electron population, say with a particle 
energy distribution harder than oc E~ 3 , given the value of 
Co, and accounting for Doppler boosting, the IC radiation 
will be mostly radiated in the GeV-TeV range, as seen by 
the observer. 

From Eq. ([ 5 ]), the total observer IC luminosity is com¬ 
puted, obtaining Tic,obs ~ 10 35 ? 7 nt erg s _1 , which is ~ 
1 r) nt% of the jet luminosity injected in the grid, and a 
~ 3ryNT% of U s h. On the other hand, neglecting relativistic 
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Fig. 5. Distribution of the ratio of fluxes of thermal energy 
to kinetic plus thermal energy. 


effects and confining the emitter to the RG wind vicinity, 
one would get Fic,obs ~ 3 x 10 32 ?7nt erg s _1 , which is ~ 300 
times smaller than the more realistic value. This shows that 
the radiation can be underestimated if the interaction re¬ 
gion is not globally taken into account. 

Figure [6] shows the structure of the observer IC emit¬ 
ter for 0 o bs l/7sh, produced in ring-like cells of volume 
27rrjArjA2j, where Arj = r max /300 and Azj = z m ax/600. 
As seen in the figure, most of the emission comes from the 
closest regions to the RG, but also from farther away, close 
to the jet shock. 

4.2. Generalization of the results 

For the adopted approximation Zj ~ bj-Rp our results cor¬ 
respond to an emitter located at relatively high heights 
(~ 10 pc) in a jet of moderate power (~ 10 44 erg s _1 ). 
These results can nevertheless be generalized through 

L [(J OC trirad/Rad ^ Lj R , (6) 

so Tic.obs will grow closer to the jet base and/or with a more 
powerful jet. In saturation, Eq. § is not valid anymore as 
^rad ^ ^nrad) and thus Eq. <© becomes 

L ic a Sbs ~ 10 2 ??NTCsh ~ 4 x 10 38 77nt erg s _1 , (7) 

i.e. it gets constant, for the same remaining parameters: yj 
and Mv w . 

The value of Tic.obs is also oc y 2 , for 0 o b s -C l/y s h and 
yj,sh 1. This relation comes from the relation 7 s h oc yj, 
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Fig. 6. Distribution of the observer IC luminosity in cgs 
units for the shocked jet region. 


strictly valid in the highly relativistic case, and the link 
between the IC luminosity in the LF (Tic.iab) and Tic,obs: 


Aic.obs/Aic.iab oc <5/ h /y 2 h « 16y s 2 h oc yf. (8) 

Therefore, the computed estimates can be analytically gen¬ 
eralized through 

Aic.obs oc y 2 J R;? g Lj,z j " 1 oc y 2 Mvy/LjZ^ 1 , (9) 

or in saturation, through 


-^IC,obs 7j 


( 10 ) 


relations that also account for different types of winds via 
Mv w . 

For a microquas ar jet, given the typical clump sizes (e.g. 
Owocki et al.||2009|), clumps with R eS ~ i?j may reach the 
jet. Also, 1C emission may take place in saturation given the 
proximity of the stellar companion, yielding luminosities as 
high as 


Lie ,obs ^ 16 y 2 h r?NTAj, 


( 11 ) 


which could be Lj for some relatively exceptional, and 
brief, events. 


5. A population of red giants within an AGN jet 

We explore in what follows whether collectively RG may 
be relevant for the non-thermal emission from AGN jets. 
About a 5 % of t he stars i n the AGN vicinity are expected 
to be RG (Syer & Ulmer 1999). In addition, there are ob¬ 
servational hints of stars interacting with AGN jets (e.g. 
Hardcastle et al. 2003 Muller et al. 20141. Therefore, the 


5 























Bosch-Ramon, V.: Non-thermal emission from standing relativistic shocks 


scenario explored in this work, and already postulated as 
the origin of high-energy emission by several authors (see 
Sect. [Tj) , could be actually a rather common phenomenon. 

We take the nearby, gamnra-ray emitting AGN M87 as 
an example. With ~ 10 11 stars i n the i nner ~ 1 kpc of the 
galaxy (see fig. 7 in Gebhardt & Tho mas|2 009), the number 
of RG within one jet should be about 


iV RG - 0.05 x 10 11 x 3a?/4 « 4 x 10 7 , 


( 12 ) 


where we take aj ~ 0.1, and 3a?/4 is the jet to total volume 
ratio. 

Assuming an RG density n RG cx 1 / , an average 

M w v w ~ 10 25 g cm s -2 , the relation from Eq. (J6J) , neg¬ 
ligible synchrotron radiation, a constant 7 j = 10, and 
Lj = 10 44 erg s _1 , each interaction would produce in 
gamma rays 


7uC,obs(2j) ~ Ai Go bs(~min) X (z m i n /Zj ), (13) 

where z m i n w8x 10 15 cm is found through Lic,obs(10 pc) = 
10 35 ?/nt erg s^ 1 , and fixing the z m i n -lunrinosity to the sat¬ 
uration value: Eic.obs^min) = ^Jcfobs- All this allows us to 
derive 

/ Zmax 

«-RG(^j)AlC,obs(2j) dvj = (14) 

min 


= M 


RG 


3 -r 

2-C, 


' 2-f _ 2-5 
^max /6 min 
3—£ _ 3-£ 
. £max Z • 


Lie .obs (-2-min ) X Z n 


where dv j is the jet differential volume. 

The dominant contribution to AjQ obs come from 
the inner-jet regions, or from galactic scales, depending 
on whether £ is > or < than 2. For £ ss 2, the z m i n , max- 
dependence of Ljg obs i s wea -k, but the external jet re¬ 
gions dominate. For this ^value^ Eq. (14) gives Ajcfobs ~ 
5 x 10 41 77 nt erg s^ 1 , implying that the contribution from 
all RG within 1 kpc could be up to a factor of a few ab ove 
the persistent GeV-TeV lumin osities in M87 (e.g. Abdo 
et al.|[2009 Aleksic et al. 2012). It seems unlikely however 
that the acceleration efficiency will be ?7 nt ~ 1, and we 
may be overestimating the average value of Mv w . On the 
other hand, we do not take into account the asymptotic gi¬ 
ant branch (AGB) and massive star contribution, both of 
which should be added to that of RG. Non-stellar objects 
such as clouds could also be relevant. Finally, the value of 
A* could be significantly higher, at least in some stars. In 
conclusion, AR?* bs may explain the persistent GeV-TeV 
emission in M87, and could be significant for AGN jets in 
general. 

Note that the viewing angle with respect to th e M87 
jet is not zero, but rather ~ 20° (e.g. Acciari et al.| [2009). 
However, given the obtained values of y s h (<; 4), an an¬ 
gle of 20° is small enough for the condition 0 o b s -C l/7sh 
to marginally apply. Here, the gamnra-ray luminosity for 
M87 has been computed using the proper viewing angle, 
although taking 0 o b s = 0 yields very similar results. 


Gebhardt & Thomas 


1 A slope £ ~ 2 can be derived from_ 

p009) for kpc scales, becoming harder closer to the centre of 


6. Discussion 


The emission coming from the jet in Blazar AGN will be 
more boosted than that from the downstream region of a 
standing shock. Thus, strong jet activity could mask the 
putative jet-obstacle emission. However, for slightly mis¬ 
aligned, relatively nearby AGN, or those without signifi¬ 
cant non-thermal activity intrinsic of the jet, the contri¬ 
bution from populations of stars interacting with the jet 
could be important, ~ 0.1 — 1% of Aj (but strongly depen¬ 
dent on ?7nt, 7j, and 0 O b s ). This collective emission should 
be constant unless £ were significantly softer than 2, which 
does not seem realistic. There could also be rare flaring 
events when a dense cloud, an AGB, a Wolf-Rayet or a lu- 


close to its base (e.g. Araudo et al.12010 

Barkov et al. 2010 

Hubbard & Blackman 2006 

Araudo et a 

. 2013). In the case 


the extragalactic case and fluctuations may be more com¬ 
mon, with the arrival of larger clumps leading to brighter 


events ( Owocki et al.||2009 Araudo et al.|[2009 ). This may 


also hap 

pen in binary systems hosting a non-accreting pul- 

sar (see 

Bosch-Ramon 

2013 

Paredes-Fortuny et al. 

2015 


and references therein). 

The obtained numerical results have been used to char¬ 
acterize the effective section of the RG wind as an obstacle 
for the jet, the distribution of the non-thermal particles 
through the constant 77 nt, and the effect of Doppler boost¬ 
ing in the shocked flow for 7 j 1. Previous research fo¬ 
cussed on a region downstream of the jet shock next to the 


obstacle and with a similar size (e.g. Araudo et al. 2009 


2010 Barkov et al. 2010 2012b), and Doppler boosting 


effects were not considered as the flow is mildly relativis¬ 
tic there. As shown here, both approximations can lead to 
the underestimation of the non-thermal emission from jet- 
obstacle interactions. 


Our work has some caveats. The simulation, in axisym- 
metry, without thermal cooling nor magnetic field, with 
modest resolution, and with a simple prescription of the 
RG wind, is not very realistic, although it captures the 
basic features of the shocked region of a hydrodynamic jet. 
The value of / ra d is a rough characterization of the radiative 
efficiency for IC, and our model of the IC emission is very 
simplified: it assumes IC only in Thomson, and accounts 
only for the electron energy at the maximum of the IC 
cross section. This simplification introduces an uncertainty 
in the predicted luminosity by a factor of a few. On the 
other hand, our radiation calculations account for the main 
features of IC scattering, and so we expect our results to be 
correct within an order of magnitude. Note that the value 
77 NT cannot be derived from first principles, and despite 
7 ?nt > 0.1 are not uncommon in relativistic jet sources, 
? 7 nt 1 are also possible. It is also difficult to constrain 
the magnetic field strength in jets, which affects the calcu¬ 
lation of synchrotron radiation. More accurate calculations 
will be presented in future work. 
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